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Abstract 



We study the behaviour of the 2d Ising model in the symmetric high temperature 
phase in presence of a small magnetic perturbation. We successfully compare the 
quantum field theory predictions for the shift in the mass spectrum of the theory 
^ ' with a set of high precision transfer matrix results. Our results rule out a prediction 

■ for the same quantity obtained some years ago with strong coupling methods. 



1 Introduction. 



Despite its apparent simplicity and the fact that since 1944 the exact expression of the free 
energy along the h = axis is known exactly the two dimensional Ising model is still 
an endless source of interesting and challenging problems P]-|24j. 

In particular in these last few years a renewed interest has been attracted by the study 
of the model in its whole complexity, i.e. in presence of both magnetic and thermal per- 
turbations. There are two main reasons behind this renewed interest. On one side recent 
progress in 2d quantum field theory (in particular, but not only, in the framework of per- 
turbed CFT's and S-matrix integrable models) allowed to obtain a host of new important 
results on correlation functions, amplitude ratios and more generally on perturbations 
around integrable directions C)n the other side the appearance of new power- 

ful algorithms (both for montecarlo simulations and for transfer matrix studies) and the 
increasing performances of computers allowed to precisely test a lot of the above predic- 
tions [in]-|211- This paper, which is devoted to the study of Ising model in presence of a 
small magnetic perturbation around the integrable thermal axis, is a further step in this 
direction. It can be considered as the natural continuation of jTHI in which, with similar 
tools, the opposite setting of a small thermal perturbation around the magnetic integrable 
line was studied and it is part of our more general ongoing project [IH]-|2l| on the study 
and characterization of the 2d Ising model in the whole thermal and magnetic plane. A 
further reason of interest for the present analysis is the theoretical estimate recently ob- 
tained by Fonseca and Zamolodchikov |3] for the first order magnetic correction to the 
mass spectrum of the model. This estimate turns out to be in sharp disagreement (more 
than one order of magnitude) with a previous estimate of the same constant obtained more 
than 20 years ago with strong coupling methods in [25j . This disagreement is very puzzling, 
since strong coupling estimates in the 2d Ising model are usually highly reliable. One of 
the goals of the present work is to fix this disagreement testing these two estimates with 
a high precision transfer matrix calculation directly in the 2d Ising model. This paper is 
organized as follows. In sect. |21we first recall some known results on the 2d Ising model 
and then concentrate on the definitions and theoretical predictions for the quantities which 
we measure in our transfer matrix analysis. Sect. |21is devoted to a short discussion of the 
transfer matrix method while in the last section we discuss our results and compare our 
findings with the theoretical predictions. 



2 The model. 



In this section we briefly recall some well known results on the Ising model. A more detailed 
description can be found in the reviews - |29p . 

^The book is the standard reference for the lattice Ising model. An updated version by one of the 
authors can be found in [27| . A recent thorough review of the field theoretic approach to the model can 
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2.1 The lattice version of the model. 



The Ising model in a magnetic field is defined by the partition function 

where the field variable takes the values {±1}; n = (no, ^i) labels the sites of a square 
lattice size Lq and Li in the two directions and lattice spacing . [n, m) denotes nearest 
neighbor sites on the lattice. In the following we shall treat asymmetrically the two direc- 
tions. We shall denote no as the compactified "time" coordinate and ni as the space one. 
The number of sites of the lattice will be denoted by = LqLi. The lattice extent in the 
transverse ("time") direction will be denoted as Lq. Following the standard notation we 
define hi = jSH thus the partition function becomes: 

Z{hl) = ^ g/^E(„,m) '^mO'^+'^iEn'^" _ (^2) 
<Ti=±l 

As it is well known the 2d Ising model has a second order phase transition at /i^ = 
and 13 = 13c given by: 

j3 = j3c = ]^\og{^/2 + I) = 0.4406868... 

In the following we shall be interested in the high temperature phase of the model (i.e. 
j3 < (3c) in which the Z2 symmetry is unbroken. 



2.2 Continuum theory 

In the continuum limit the model is described by the action: 

A = Ao — T / <f'xe{x) — h I cfxa{x), (3) 



where cr{x) and e{x) denote the magnetic and thermal perturbing operators respectively 
while h and r are dimensional couplings measuring the magnetic field and the deviation 
from critical temperature. A standard scaling analysis shows that the scaling behaviour of 
these two couplings is 



r ~ = m . 



h ~ m^-^- = m^^/^ (4) 

where m ~ 1/^ (^ being the correlation length of the model) denotes the mass scale 
associated to the breaking of scale invariance away from criticality, = 1 and x^ = 1/8 

be found in 

^ Since the lattice spacing will play no role in the following we shall set a = 1 in the rest of the paper. 
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are the scaling dimensions of the energy and spin operators respectively and denotes 
the conformal invariant action at the critical point. In the Ising case, which is known to 
be described by a free massless Majorana fermion, this can be written explicitly as: 




(5) 



where d = dz = {di — 182) 12 and d = = {di + 182) and the operators ifj and are the 
two components of the neutral Majorana fermion. 

The thermal perturbing operator (which coincides with the energy operator) can be 
written as 



which allows to identify it as the mass term of the free fermionic action, in agreement with 
the first of the scaling equations (jH). 

By combining the two scaling equations Q one can see that the field theory Q describes 
a one-parameter family of renormalization group trajectories flowing out of the critical 
point at r = /i = and labeled by the dimensionless quantity^ 



The Ising field theory can be solved exactly in the two limiting cases of /i = and 
r = 0. In the first case, as we have seen above, it is simply the theory of free massive 
fermion, the mass being proportional to |r|. In the second case (which corresponds in the 
lattice discretization to the magnetic perturbation with (3 fixed to the critical value j3 = f3c) 
A. Zamolodchikov was able to show that Q with r = is a complicated but integrable 
quantum field theory j2] of eight interacting particles. 

For generic values of h and r exact integrability is lost but notwithstanding this some 
important theoretical results can all the same be obtained. Let us see a few of them. 

• Equation of state. 

With a combination of transfer matrix techniques and analytic continuation of suit- 
able parametric representations it is possible to construct a very precise expansion 
for Helmholtz free energy of the the model and from it of the equation of state in the 
whole critical region in the (r, h) plane [21]. From this expression precise predictions 
for several universal amplitude ratios can be obtained (see ^T] for details). 

• The spectrum of the model. This important result was conjectured for the first 
time in a seminal paper of Mc Coy and Wu jHO] and recently re-understood in a field 

■^To avoid confusion let us stress that we follow the definition for 77 of jSH] which is different from the 
one adopted in |4|. 



(6) 



r 



(7) 



V 



/i 1 8/15 ■ 
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theoretic language (see [21] for a thorough discussion). The scenario which emerges 
is the foUowing. For t] = — oo (i.e. in the low temperature broken symmetry phase) 
the spectrum starts with a cut. In this limit the natural degrees of freedom of the 
model are kinks interpolating between the two degenerate vacua which are non-local 
with respect to the spin degrees of freedom of the model. As the magnetic field is 
switched on (i.e. for large and negative values of rj) this cut breaks up in a host of 
particles which may be understood as bound states of the above mentioned kinks. As 
T] is increased from —oo to finite negative values the number of stable particles (i.e. 
below the pair creation threshold rUn < 2mi) decreases. Exactly at r] = only three 
such particles are left (but, only for this value of rj, due to the exact integrability of 
the model five more particles turn out to be stable even if they are above threshold). 
The number of stable particles continues to decrease as rj increases, until a single 
particle is left at large enough 77. There is a wide region of values of rj before the 
positive thermal axis (i.e. 77 = +00) is reached in which the spectrum contains only 
one stable particle. This is the region in which we shall perform our analysis in the 
following. 



00 




Figure 1: Ising Model's masses, the value of the masses has been obtained in [3] 



• Small perturbations around the integrable lines. The most interesting result 
for our present purposes is that in the nearby of the two integrable lines a few exact 
results can be obtained using perturbative methods. In particular one can evaluate 
the corrections to the masses and to the ground state energy due to the perturbing 
operator [7]. Again we have two possible situations: 

a] A small thermal perturbation of the Zamolodchikov's integrable model along 
the r = magnetic axis. Estimates for the first order mass correction where 



4 



obtained in [7j and successfully tested both on the lattice, using transfer matrix 
methods JH] and directly in the Ising field theory by a numerical diagonalization 
of the Hamiltonian on a conformal basis of states [Zj. 

b] A small magnetic perturbation of the free fermionic theory along the h = axis. 
This case was recently studied by Fonseca and Zamolodchikov j3] which were 
able to obtain, using the Ward identities of the model, the matrix elements of 
the product cr(x)o"(x') between any particle states. This is exactly the ingredient 
which is needed to obtain the sought for mass correction. The aim of our paper 
is to test this result with a transfer matrix analysis directly on the lattice model. 
We shall devote the next subsection to a detailed discussion of this term. 

2.3 The mass correction 

Following the above discussion in the r > phase, for small value of h there is a single 
stable particle in the spectrum. Let us call m(r, h) its mass. For symmetry reasons the 
correction to the mass of the particle when the magnetic field is switched on at r > must 
be proportional to h"^. Scaling arguments (see the discussion in the previous section) then 
fix the dependence on r to be r~^^/^. Thus we expect the following behaviour: 

m(r, h) = m(r, 0)(1 + a/iV"^^/^) (8) 

From a numerical point of view it is useful to rewrite this dependence as a function of the 
unperturbed mass m(r, 0): 

m(r, h) = m(r, 0)(1 + a/i^m(r, 0)"^^/^) (9) 

Recently, in a remarkable paper Fonseca and Zamolodchikov |lj were able to evaluate 
analytically this constant which turns out to be 

a = 10.7619899(1) (10) 

The above correction can also be rewritten in a form which will be useful in the following: 

m2(r, h) = m^T, 0) + fi/iV"^/^ (11) 

These three constants, which are actually the same constant written in different units, can 
be easily related among them. In fact, thanks to the fact that the Ising quantum field 
theory is solved exactly along the thermal line, the relation between r and m(r, 0) is know 
exactly^: m = 2txt. This allows to write the following relations: 

~a = 2a(27r)-^/^ (12) 
a = a(27r)-i5/^ (13) 

"'This relation, as all the ones which follow in this section actually depends on the conventions which 
one chooses for the structure constants. In this paper we follow the so called "conformal normalization" 
i.e. Ci, = Cl, = 1 . 
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2.3.1 Theoretical estimates for mass correction 

Up to a few months ago the only existing estimate for the constant a was the one reported 
in [Snj, obtained performing the continuum limit extrapolation of the strong coupling 
expansion of the perturbed two point correlator ^: 

1 1 

a = (- + -)(2V2)3 = 0.73700676... (14) 

4 TT 

This estimate was never tested with numerical simulations or transfer matrix methods. 
Further, the impressive discrepancy between ()10p and the above estimate was one of the 
major reasons which prompted us to perform the present transfer matrix analysis. 

2.4 Correction in the free energy. 

Similarly to the mass case one can define the deviation from the unperturbed value of 
the free energy JF (which is nothing else than the lowest eigenvalue of the transfer matrix 
spectrum). In this case it is easy to obtain a theoretical estimate of this deviation by 
noticing that for symmetry reasons the first nonzero correction must be again at order h? 
and is exactly given by 1/2 times the magnetic susceptibility 

^(r, h) = T{t, 0) + yx+(r, 0) + 0{h') ^ T{t, 0) + ^h\''l^ + ... (15) 

Where in the last term of the above equation r+ denotes the susceptibility amplitude and 
we have inserted the known dependence on r of the susceptibility. An important role in the 
following analysis is played by the fact that the value of r_|_ in the the continuum theory 
is known exactly (see for instance p9|): 

r+ = 0.148001214... . (16) 

2.5 Universal amplitude ratios. 

A major problem in comparing the above values, obtained in the continuum limit QFT, 
with our numerical estimates in the lattice model is that we must convert the continuum 
limit parameters r and h into their lattice counterparts t and hi. A nice way to avoid 
this problem is to combine the above constants in a suitable universal ratio which is thus 
independent from the particular realization of the underlying quantum field theory. There 
are two natural choices: 

^We have normalized the result of |^ in units of the continuum field theory discussed above, thus it 
may be directly compared with the one of eq. ((Tn|l . 
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Following j29j we shall study the combination 

6 

Jj— >oo SJ' 



= liin ^ (17) 



where 

6m^ = m^{T,h) -m^{T,0) (18) 
6T = J^{T,h)-J^{T,0) (19) 



A direct substitution gives: 



A second option is: 



where 



Smf 4a 
^^^ = f7(2^ (20) 



_ 6m M^{-r) 
' m{T,0)M^{T,h) ^ ' 



6m = m{T,h) -m{T,0) (22) 



while M{T,h) denotes the magnetization and Ms(— r) the spontaneous magnetiza- 
tion. A direct substitution of the various terms leads to 



R, = a{27r)-''/'^ (23) 

where B denotes the amplitude of the spontaneous magnetization whose value is 
known exactly in the continuum limit theory (see for instance |29j) 



B = 1.70852190. (24) 

The expected values for these ratios, using the Fonseca-Zamolodchikov estimate for the 
constant a are 

Ri = 11.66467... R2 = 1.4568962... (25) 

while with the strong coupling estimate for a we should expect: 

Ri = 0.7988246... R2 = 0.09977173... (26) 



2.6 Lattice amplitudes. 

This is all we can do for a generic realization of the Ising quantum field theory. However 
in the particular case of the 2d square lattice realization thanks to the fact that the model 
is exactly solved also on the lattice for /i = we know the value of several of the above 
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amplitudes directly in the lattice realization. The mass of the theory in the transfer matrix 
geometry is 

rrii = - log V ^'^^^ 

with V = tanh/5. The index i recalls that this is the lattice estimate of the mass of the 
model. Eq. (j27j) leads to the following expansion in the vicinity of the critical point 

m,(t,0) ~4/?,t (28) 



where t is the reduced temperature defined as 

1. 



(=^. (29) 



The spontaneous magnetization is given by 



Ms = (l- ' (30) 

which gives: 

M, ~ (-8V2/3et)i (31) 

For the magnetic susceptibility x+(^) there is no exact expression. However this ob- 
servable has been studied extensively in the past years, firstly by McCoy et al. PI] and 
more recently by Orrick et al. 32] using strong coupling expansions (notice that the first 
term in the expansion in powers of t can be evaluated exactly). One finds: 

x+(t) = Co,+r^/^ + Ci,+r3/^ + ... (32) 

where 

Co,+ = 0.962581732... 

Ci^^ = 0.074988153... . (33) 

Plugging these amplitudes (together with the lattice estimate of the constant a which we 
shall discuss in the following section) it is possible to obtain the two ratios Ri and R2 
discussed in the previous section. 

Another possible use of these results is to construct the explicit relation between t and 
r and similarly between he and h. This will allow us to test separately the two predictions 
eq.® and (fT3j). Comparing the lattice and continuum values for the mass we find: 

r = i^5(i±^i ^ i, (34) 

TT 
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while comparing the susceptibility amplitudes 

Co,+ = K'/'kl r+ (35) 

we find 

h = 25/^8ei/«A-3/2/,^ = h, (36) 
(where A = 1.28242712 ... is the Glaisher constant). 

With these results at hand we may write an explicit prediction for the constant a 
directly on the lattice. Defining the lattice version of a as: 

ae = klk;^^'^ a (37) 

we find using the value for a obtained by Fonseca and Zamolodchikov |1] 

= 7.56977... (38) 

and, using the strong coupling one discussed in pK] : 

= 0.518396... . (39) 

In the following we shall use this last result and shall directly extract from the transfer 
matrix data the value of which we shall then combine with the other amplitudes in order 
to obtain the two universal ratios -Ri and i?2- 



3 The transfer matrix analysis 

For the analysis of the transfer matrix data we followed the same procedure of our previous 
works IT7| 1^ which can be schematically sketched as follow. 

• We evaluate the first few eigenvalues of the transfer matrix for the values of tem- 
perature and magnetic coupling in which we are interested for lattice widths in the 
transverse direction in the range Lg{9 — 21}. We construct from them the observ- 
ables in which we are intested (mass and free energy) and extrapolate their infinite 
lattice limits using the recursive procedure discussed in detail the appendix of fFI\ . 

• For each values of t and we construct the combinations 5m and and 5JF 
defined in eq.(I221), (HHI) and 

• We fit 5m? , 5m and 5J^ assuming the known scaling behaviour for these quantities. 
As we shall see the first few terms of the scaling functions will be enough to obtain 
stable results for the fits within the precision of our data. 

In table (Q) we report the parameters and the various settings we used in our simulations, 
while in tables ()2|3j) we report a sample of our data to allow the interested reader to 
reproduce our analysis. 
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t 


hi 


_^ 




n n 

U.U 


Q 


0.350 


0.0003 


10 


0.355 


0.0005 


11 


0.360 


0.0008 




0.365 


0.0010 




0.370 


0.0013 






0.0015 


19 




0.0018 


20 




0.0020 


21 



Table 1: Parameters value. 



L 


l/m{t, hi) 


T{t,hi) 


Mt{t, hi) 


15 
16 
17 
18 
19 
20 
21 


2.565917754587 

2.565918168084 

2.5659182662035 

2.565918247605 

2.565918227627 

2.565918215050 

2.565918209646 


0.830031898983 
0.830031896272 
0.830031895299 
0.830031894746 
0.830031894142 
0.830031894612 
0.830031893614 


0.020167928578 
0.020167777385 
0.020167731420 
0.020167716025 
0.020167710600 
0.020167708951 
0.020167708217 


oo 


2.56591821(1) 


0.830031894(1) 


0.020167708(1) 



Table 2: Example of data used for the extrapolation {hi = 0.0013, (3 = 0.350). 

4 Discussion of our results. 

Now we are in the position to compare the lattice results, coming from the transfer matrix 
approach, with the field theoretic estimates of the ratios Ri, R2. To this aim, the first step 
is to extract the leading behaviour of 6m, dm?, and Mi{t, hi) from the data obtained 
by means of the transfer matrix. 

As a preliminary check of our procedure we estimated from the transfer matrix data 
the mass at zero magnetic field which in this geometry is known exactly and is given by 
eq.(|77|). We report in tab 01 the comparison which turns out to be fully satisfactory. 

Then we concentrated in the evaluation of brr? . Symmetry considerations suggest the 
following expansion in even powers of hi for this quantity. 

6m' = h{t)h'i + f2{t)ht + . . . (40) 

It turns out that, within the precision of our data and given our choice of values of hi, it 
is enough to take into account only the first two terms in this fit, the term proportional to 
h^ being completely negligible. For each fixed value of t we extract from the fit the best 
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fit values for and f2{t)- Since we are interested in tlie leading behaviour of Sm"^, in 
the following we shall concentrate in the study of the function fi{t). The functional form 
of fi{t) can be easily constructed using the same RG and CFT arguments developed in 

m-m. m-- 

where ae is the amplitude in which we are interested. 
Following the same fitting procedure used in we obtain: 

a£ = 5.61 ±0.01 (42) 

which gives 

ae = 7.56 ± 0.01 (43) 

in complete agreement with the Fonseca-Zamolodchikov estimate of the amplitude a (see 
eq. dSHD). 



The amplitude of the leading behaviour of 5m can be immediately inferred from that 

of 

= ht-^^'%l + ... (45) 

plugging in the last definition the above obtained value of an we obtain 

= 1.591 ± 0.002. (46) 

As a consistency check, we computed the same amplitude directly from the fit of the 
transfer matrix data. The final result is 

he = 1.592 ± 0.003 (47) 

which is completely equivalent to the previous estimate. 

To estimate the vacuum energy on the lattice let us notice first of all that 5jF can be 
expanded as 

5J'=^X+it)h',+g2{t)ht + ... (48) 

where the function x+{t) is the magnetic susceptibility of the model in the high-temperature 
phase (3 < (3c and h = already discussed in section 12.61 Hence the exact estimate of the 
amplitude of the free energy on the lattice is given by 

Q = ^ = 0.48129086 .... (49) 
11 



hi 


/3 = 0.345 


/3 = 0.350 


/? = 0.355 


/3 = 0.360 





2.42236580(1) 


2.56721580(1) 


2.72886532(1) 


2.91043999(6) 


0.0003 


2.42231285(1) 


2.56714664(1) 


2.72877368(1) 


2.91031661(4) 


0.0005 


2.42221873(1) 


2.56702370(1) 


2.72861080(1) 


2.91009730(5) 


0.0008 


2.42198936(1) 


2.56672412(1) 


2.72821391(1) 


2.90956300(5) 


0.0010 


2.42177770(1) 


2.56644770(1) 


2.72784773(1) 


2.90907012(4) 


0.0013 


2.42137222(1.5) 


2.56591821(1) 


2.72714641(1) 


2.90812629(4) 


0.0015 


2.42104331(1) 


2.56548878(1) 


2.72657770(2) 


2.90736112(3) 


0.0018 


2.42046226(1) 


2.56473025(1) 


2.72557340(1) 


2.90601025(2) 


0.0020 


2.42001655(1) 


2.56414851(1) 


2.72480337(1) 


2.90497486(3) 



(a) Inverse mass l/m{t, hi) 


he 


(3 = 0.345 


/? = 0.350 


(3 = 0.355 


/3 = 0.360 


0. 

0.0003 
0.0005 
0.0008 
0.0010 
0.0013 
0.0015 
0.0018 
0.0020 


0.825668076(1) 
0.825668712(1) 
0.825669843(1) 
0.825672600(1) 
0.825675145(1) 
0.825680021(1) 
0.825683978(1) 
0.825690973(1) 
0.825696341(1) 


0.830018782(1) 
0.830019481(1) 
0.830020722(1) 
0.830023748(1) 
0.830026541(1) 
0.830031894(1) 
0.830036238(1) 
0.830043915(1) 
0.830049807(1) 


0.8344669273(3) 

0.834467698(1) 

0.834469068(1) 

0.834472408(1) 

0.834475490(1) 

0.834481397(1) 

0.834486189(1) 

0.834494660(1) 

0.834501162(1) 


0.839014965(1) 
0.839015821(1) 
0.839017342(1) 
0.839021049(1) 
0.839024471(1) 
0.839031028(1) 
0.839036348(1) 
0.839045750(1) 
0.839052965(1) 


(b) Free energy J^{t, hg) 


he 


13 = 0.345 


(3 = 0.350 


13 = 0.355 


13 = 0.360 


0.0003 
0.0005 
0.0008 
0.0010 
0.0013 
0.0015 
0.0018 
0.0020 


0.004241449(1) 
0.007068848(1) 
0.011309247(1) 
0.014135509(1) 
0.018373549(1) 
0.021197803(1) 
0.025432177(1) 
0.028253554(1) 


0.004655977(1) 

0.00775956(1) 

0.012414204(1) 

0.015516337(1) 

0.020167708(1) 

0.023267128(1) 

0.027913548(1) 

0.031009083(1) 


0.005138273(1) 

0.00856335(1) 

0.01369968(1) 

0.017122651(1) 

0.022254596(1) 

0.0256738415(5) 

0.030798990(1) 

0.034212903(1) 


0.00570442(1) 

0.00950676(1) 

0.01520845(1) 

0.019007841(5) 

0.024703401(2) 

0.028497574(2) 

0.034183633(2) 

0.037970351(1) 


(c) Magnetization M{t, hi) 



Table 3: Extrapolated values of different operators as function of the magnetic coupling 
and the temperature. 



We also compared the previous exact results (quoted in sect. 12 .6^ with our numerical 
estimates. We found 

= 0.962 ±0.002 
= 0.075 ±0.001 

in complete agreement with ()H3j) . 

Finally, we need the amplitude di of the scaling behaviour of the magnetization. It can 
be easily computed (exactly) from the definition of M^it, hi) 

Mi{t, hi) = — = Co,+ + ... (50) 
orii 
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Exact mass formula 


Transfer matrix 


0.345 


2.4223658021832 


2.42236580(1) 


0.350 


2.5672157958004 


2.56721580(1) 


0.355 


2.7288653104825 


2.72886532(1) 


0.360 


2.9104399518676 


2.91043999(5) 


0.365 


3.1158911582069 


3.115891(1) 


0.370 


3.3502882801313 


3.350287(1) 



Table 4: Comparison of our TM extrapolation with the exact result for the invers mass at 
zero magnetic field 



leading to 

= Co,^ = 0.962581732 .... (51) 

Since the transfer matrix approach allows also the numerical computation of the magneti- 
zation, we can extract the numerical estimate of di 

de = 0.961 ± 0.004 (52) 

which, once again, agrees with the previous results. We also stress that the zero-field 
quantities and Mg involved in the ratio R2 are given by their exact expressions quoted 
in section 

The final step is then to compute the lattice estimates of the ratios Ri, R2 

R^ = ^ = 11.66 ±0.02 
Q 

i?2 = 2-9/Vc"^/^-^ = 1.456 ±0.002 (53) 
dj 

which agree with the field theoretic predictions obtained using the Fonseca and Zamolod- 
chikov estimate of the constant a (see eq. fE3|) ). 

The precision of our numerical computations allows to rule out unambiguosly the strong 
coupling estimates of these universal ratios. 

It would be interesting to investigate the reasons of this failure, which is even more 
surprising given the usual high reliability of strong coupling results in two dimensional 
statistical models (and in particular in the Ising model). A possible explanation is that 
the result of |25j is based on a peculiar mutipole expansion of the perturbed two point 
function, which was proved to have good converegence properties when tested in the case 
of the H = susceptibility (see [21]) but most probably has a less good behaviour if one 
is interested in the mass shift. 



13 



Acknowledgments. We would like to thank G. Delfino for useful discussions. This 
work was partially supported by the European Commission TMR programme HPRN-CT- 
2002-00325 (EUCLID). The work of P.G. is supported by the COFIN "Teoria dei Campi, 
Meccanica Statistica e Sistemi Elettronici" . 



References 

[1] L. Onsager, Phys. Rev. 65 (1944) 117. 

[2] A. B. Zamolodchikov, Adv. Stud. Pure Math. 19 (1989) 641, 
A. B. Zamolodchikov, Int. J. Mod. Phys. A 4 (1989) 4235. 

[3] P. Fonseca and A. Zamolodchikov, J. Stat. Phys. 110 (2003) 527, 
iarXiv:hep-th/0112167j . 

[4] P. Fonseca and A. Zamolodchikov, |arXiv:hep-th/ 0309228 1 . 

[5] V. A. Fateev, Phys. Lett. B 324 (1994) 45. 

[6] V. Fateev, S. Lukyanov, A. B. Zamolodchikov and A. B. Zamolodchikov, Nucl. Phys. 
B 516 (1998) 652. arXiv:hep-th/9709034, . 

[7] G. Delfino, G. Mussardo and P. Simonetti, Nucl. Phys. B 473 (1996) 469 
|arXiv:h ep-th/ 9603011i . 

G. Delfino and G. Mussardo, Nucl. Phys. B 455 (1995) 724. |arXiv:hep-th/9507010| . 

[9] G. Delfino, P hys. Lett. B 419 (1998) 291 [Erratum-ibid. B 518 (1998) 330]. 
|arXiv:hep-th/9710019|. 



[10] G. Delfino, J. Phys. A 34 (2001) L161. |arXiv:hep-th/0101180 



[11] A. Leclair and G. Mussardo, Nucl. Phys. B 552 (1999) 624. |arXiv:hep- th/9902075l 

[12] G. Mussardo, J. Phys. A 34, 7399 (2001). |ar Xiv:hep-th/0103214| . 

[13] P. Grinza and N. Magnoli, J. Phys. A 36 (2003) L509 |arXiv:hep-th/0306100|. 



[14] P. Calabres e, M. Caselle, A. Ceh, A. Pelissetto and E. Vicari, J. Phys. A 33 (2000) 
8155 arXiv:hep-th/0005254 . 



[15] M. Caselle, P. Grinza, R. Guida and N. Magnoli, J. Phys. A 37 (2004) L47 
|arXiv:h ep-th/03l36086 1 . 

[16] P. Grinza and A. Rago, Nucl. Phys. B 651 (2003) 387 ar Xiv:hep-th/0208016| . 

[17] M. Caselle and M. Hasenbusch, Nucl. Phys. B 579 (2000) 667. |arXiv:hep-th/9911216| . 



14 



M. Caselle, P. Grinza and N. Magnoli, Nucl. Phys. B 579 (2000) 635. 
|arXiv:hep-th/9909065|. 

M. Caselle, P. Grinza and N. Magnoli, J. Phys. A 34 (2001) 8733. 
|iarX iv:hep-tli/ 0103263|. 

M. Caselle, M. Hasenbusch, A. Pelissetto and E. Vicari, J. Phys. A 35 (2002) 4861 
arXiv:cond-mat /0106372 1 . 

M. Caselle, M. Hasenbusch, A. Pelissetto and E. Vicari, J. Phys. A 34 (2001) 2923 
arXiv:cond-mat / 001 1305 1 . 

M. Caselle, M. Hasenbusch, A. Pelissetto and E. Vicari, J. Phys. A 33 (2000) 8171 
^arXiv:hep-th/0 003049 1 . 

M. Caselle and M. Hasenbusch, Nucl. Phys. B 639 (2002) 549 |arXiv:hep-th/ 0204088 . 
P. Grinza and A. Rago, J. Phys. A 37 (2004) L4701, |arXiv:cond-mat7"0210046] . 
S. S. Salihoglu, Phys. Rev. D 24 (1981) 3255. 

B. M. McCoy and T. T. Wu, The two dimensional Ising Model, (Harvard Univ. Press, 
Cambridge , 1 9 73 ) ; 

B. M. McCoy, in Statistical Mechanics and Field Theory, eds. V.V. Bazhanov and C.J. 
Burden, (World Scientific, Singapore, 1995). 

G. Delfino, Phys . Lett. B 419 (1998) 291 [Erratum-ibid. B 518 (2001) 330] 
||arXiv:hep-th/9710019|. 

G. Delfino, J. Phys. A 37 (2004) R45, |arXiv:hep-th/0312119 



B.M. McCoy and T.T. Wu, Phys. Rev. D 18 (1978) 1259. 

T. T. Wu, B. M. McCoy, C. A. Tracy and E. Barouch, Phys. Rev. B 13 (1976) 316. 

W. P. Orrick, B. Nickel, A. J. Guttmann, J. H. H. Perk, J. Statist. Phys. 102 (2001) 
795-841; and W. P. Orrick, B. G. Nickel, A. J. Guttmann, J. H. H. Perk Phys. Rev. 
Lett. 86 (2001) 4120-4123. 



15 



